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Abstract 

We compute a parametric description of the totally mixed Nash equilibria of a 
generic game in normal form with pre-fixed structure. Using this representation, 
we show conditions under which a game has the maximum possible number of this 
kind of equilibria. Then, we present a symbolic procedure that allows us to describe 
and estimate the number of isolated totally mixed Nash equilibria of an arbitrary 
game. Under certain assumptions, the algorithm computes the exact number of these 
equilibria. 

1 Introduction 

Noncooperative game theory is used to model and analyze strategic interaction situations. 
Among its most outstanding applications, we can mention the fundamental role this theory 
has played in economics (see, for example, the classical reference book |31J). Moreover, 
game theory has also been applied to politics, sociology and psychology, and to biology 
and evolution as well. 

One of the main concepts in this theory is a Nash equilibrium, which consists in a 
situation in which no player can increase his payoff by unilaterally changing his strategy. 
Since within this theory the players cannot communicate in order to decide a simultaneous 
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change of strategies, in a Nash equilibrium the game stabilizes. In [22], it is proved that 
any noncooperative game in normal form has at least one Nash equilibrium. However, the 
proof is not constructive and does not give any information about the existence of more 
than one Nash equilibrium. The question posed is how to compute algorithmically Nash 
equilibria and to determine the number of them in a given game. 

Nash equilibria of noncooperative games in normal form can be regarded as real solu- 
tions to systems of polynomial equations and inequalities (see, for instance, [29, Chapter 
6]). In the case of two players, each equilibrium is the solution of a linear system of equa- 
tions, and therefore, equilibria may be found exactly by using simplex type algorithms 
(see, for instance, |15j): however, there is no polynomial time algorithm solving the prob- 
lem (see |32j). In the general case of a game with more than two players, the polynomials 
appearing are multilinear. To solve the problem of finding one equilibrium, some nu- 
merical methods have been applied successfully (for example, some methods derived from 
Scarf's algorithm, [26]). Nevertheless, sometimes it is not sufficient to compute only one 
equilibrium because, depending on the problem to be solved, not all the equilibria of a 
game are equally interesting and the methods developed to compute only one equilibrium 
do not allow us to decide beforehand whether it fulfills some additional properties or to 
compare different equilibria. 

A comparative study of different known methods for the computation of all the Nash 
equilibria of a game may be found in [7J. In [12], a new algorithm solving this problem for 
generic games by means of homotopy methods is presented, but no complexity bounds are 
shown. In addition, the application of symbolic algorithms solving systems of equations 
and inequalities over the real numbers is being studied in this context, motivated by the 
characterization of the set of all the Nash equilibria of a game as a semi-algebraic set (an 
example of this fact is the application of quantifier elimination algorithms over the real 
numbers to compute approximated equilibria in [16]; see also the survey p2]). However, 
up to now, no significative result had been obtained concerning the adaptation of these 
algorithms in order to profit from the particular properties of the algebraic systems arising 
in game theory. 

In this paper, we study totally mixed Nash equilibria, that is to say, Nash equilibria in 
which every player allocates a positive probability to each of his available strategies. Note 
that a procedure to compute these equilibria can be used as a subroutine to compute all 
Nash equilibria of the game by recursing over all possible subsets of used strategies. We 
present a symbolic method to find a parametric description of the set of totally mixed Nash 
equilibria of a generic game with a pre-fixed structure. This method is based on a symbolic 
procedure for the computation of multihomogeneous resultants with complexity polynomial 
in the degree and the number of variables of the resultant ([13]). Using the description 
previously obtained, we show conditions (given by polynomial inequalities on the payoff 
values) under which a game with the given pre-fixed structure will have the maximum 
number of such equilibria. The next step is to give a similar parametric description 
for particular games. First, we solve this problem under some genericity assumptions 
implying, in particular, that the number of totally mixed Nash equilibria of the game is 
finite, and we show how to compute this number. Then, we consider the general case, in 
which we give a parametric description of a finite set of points including all the isolated (in 
the complex space) totally mixed Nash equilibria of the game, which enables us to bound 
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the number of these isolated equilibria. All our algorithms have a complexity polynomial 
in the number of pure strategies of each player, and the number of totally mixed Nash 
equilibria of a generic game with the given structure. 
This paper is organized in the following way: 

In Section [21 we introduce some basic notions on game theory, the description of the 
algorithmic model we are going to use and a mathematical formulation for the totally 
mixed Nash equilibria of a game in normal form as the set of solutions of a system of 
multihomogeneous polynomial equations. In Section [3j the algorithmic description of the 
totally mixed Nash equilibria of a generic game is given, and it is used to find conditions 
under which a game has the maximum number of these equilibria. In Section [4] we deal 
with the totally mixed Nash equilibria of particular games. Finally, the last section is 
devoted to proving some results about the algorithmic computation of multihomogeneous 
resultants and upper bounds for their degrees which are used throughout the paper. 

2 Preliminaries 

2.1 Definitions and Notation 

2.1.1 Basic facts 

Throughout this paper Q denotes the field of rational numbers, N denotes the set of 
positive integers and No :=NU {0}. 

If if is a field, we denote an algebraic closure of K by K. The ring of polynomials in the 
variables x±, . . . , x n with coefficients in K is denoted by K[xi, . . . , x n ]. For a polynomial 
/ G K[x\, . . . , x n ] we write deg / to refer to the total degree of / and deg^. / to refer to 
the degree of / in the variable Xj. 

For n G N and an algebraically closed field k, we denote by A n (k) and ¥ n (k) (or simply 
by A™ or P n if the base field is clear from the context) the re-dimensional affine space and 
projective space over k respectively, equipped with their Zariski topologies. 

We adopt the usual notions of dimension and degree of an algebraic variety V, which 
will be denoted by dim V and deg V respectively. See for instance [27] and [9] for the 
definitions of these notions. 

2.1.2 Game theory 

In this section we present some basic game theory concepts. For a more detailed account 
on the subject we refer the reader to any standard game theory text such as [24] . 

We consider non-cooperative games in normal form; that is to say, games in which there 
is only one time step at which all the players move simultaneously without communicating 
among themselves. We will assume that there are r players in the game having n\ + 
1, . . . , n r + 1 distinct available pure strategies respectively (rti, . . . , n r G N). 

For i = 1, . . . , r: 

• cW := (c^ j r )o<j k <n k is the given payoff matrix of player i, where ^ is the 
payoff to player i if, for every 1 < k < r, player k chooses the strategy jV 
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• Xi := (xio,Xn, . . . ,Xi n .) is a vector representing a mixed strategy of the ith player, 
which is a probability distribution on his set of pure strategies (that is to say, Xy is 
the probability that the ith player allocates to his jth pure strategy). 

With these notations, for every 1 < i < r, the payoff to player i if the mixed strategies 
X\, . . . , X r are played is: 



A Nash equilibrium is a vector of mixed strategies such that no player can increase 
his payoff by changing unilaterally to another strategy while the other players keep their 
strategies fixed; that is to say, a vector of mixed strategies X\, . . . , X r satisfying 



7Tj(Xi, . . . Xi, X i+ i, . . . ,X r ) > iTi(Xi, . . . Xi,X i+ i, . . . ,X r ) VI < i < r 



for every mixed strategy X[. A totally mixed Nash equilibrium is a Nash equilibrium 
in which each pure strategy is assigned a positive probability, that is, one that satisfies 
Xij > for every 1 < i < r, < j < rii. 

2.2 Data structure, algorithms and complexity 

The algorithms we consider in this paper are described by arithmetic networks over the 
base field Q (see [33]). An arithmetic network is represented by means of a directed 
acyclic graph. The external nodes of the graph correspond to the input and output of the 
algorithm. Each of the internal nodes of the graph is associated with either an arithmetic 
operation in Q or a comparison between two elements in Q followed by a selection of 
another node. These are the only operations allowed in our algorithms. 

We assume that the cost of each operation in the algorithm is 1 and so, we define the 
complexity of the algorithm as the number of internal nodes of its associated graph. 

The objects our algorithms deal with are polynomials with coefficients in Q. We 
represent each of them by means of one of the following data structures: 

• Dense form, that is, as the array of all its coefficients (including zeroes) in a pre-fixed 
order of all monomials of degree at most d, where d is an upper bound for the degree 
of the polynomial. The size of this representation equals the number of monomials 
of degree at most d. 

• Sparse encoding, that is, as an array of the coefficients corresponding to monomials 
in a fixed set, provided that we know in advance that the coefficient of any other 
monomial of the polynomial must be zero. The size in this case is the cardinal of 
the fixed set of monomials. 

• Straight-line programs, which are arithmetic circuits (i.e. networks without branches). 
Roughly speaking, a straight-line program (or sip, for short) over Q encoding a poly- 
nomial / E Q[xi, . . . , x n ] is a program which enables us to evaluate the polynomial 
/ at any given point in Q n . Each of the instructions in this program is an addition, 
subtraction or multiplication between two pre-calculated elements in Q[x\, . . . ,x n ], 




0<ji<ni 



0<j r <n r 
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or an addition or multiplication by a scalar. The number of instructions in the pro- 
gram is called the length of the sip. For a precise definition of straight-line programs 
we refer to Definition 4.2] (see also 

2.3 Polynomial equations for totally mixed Nash equilibria 

Now, we show that the totally mixed Nash equilibria of an r-person game in normal form 
can be regarded as real solutions to a polynomial equation system (see [2£l Sec. 6.3]). 

Let us observe that, if the payoff matrix of player i is (1 < i < r), given mixed 
strategies Xi > . . . , X r , the payoff to player i is 

ir i (X 1 ,...,X r ) = ^ cf^ jr x lh . . . x rjr 

31,— Jr 

= c h-jr Xl n ■ 

0<ji <"i jl,— Ji-l ,3i+U ■ -,3r 

= ^2 Xi 0i ^it^lj • • • !-^»-l> e jl ) ^i+li 
0<ji<ni 

where is the (J + l)th vector of the canonical basis of M n ' 

Now, if (X\, . . . ,X r ) is a totally mixed Nash equilibrium, we have that > for 
every i, j. Then, in order that ir^Xi, . . . , • , -Xj+i, • • • , X n ) attains a maximum at 
Xi = (xio, . . . , xi ni ) (among all the different probability distributions), a necessary and 
sufficient condition is that 

7Tj(Xi, . . . , Xi-i, Xi + i, . . . , X r ) = TTi(Xi, . . . , Xi-i,e$- ,Xi + i, . . . , X r ) V 1 < ji < ni. 

(1) 

In fact, assuming that this condition does not hold for some 1 < i < r, let j max and 
Jmin denote the indices of pure strategies leading to maximum and minimum payoffs to 
player i respectively. Then, taking e > with e < min{xj Jmin , 1 — :% max } and changing 
to the (totally) mixed strategy X- := Xi + &(e^ ax — ) the ith player can increase his 
payoff by 

e(7Tj(Xi, . . . , e^ ax , X i+ i, . . . , X r ) - 7Tj(Xi, . . . , e^. n , . . . , X r )) > 0, 

contradicting the fact that the vector is a totally mixed Nash equilibrium for the game. 

Conversely, for a vector of totally mixed strategies (X\, . . . , X r ) such that all the identi- 
ties (HJ) hold for every 1 < i < r, we have that 7Tj(Xi, . . . , . . . , X r ) = m(Xi, . . . , Xi, . . . , X r ) 
for every mixed strategy X[ and so, player z (1 < i < r) cannot increase his payoff by 
unilaterally changing his strategy, which implies that (X±, . . . ,X r ) is a Nash equilibrium 
for the game. 

We conclude that the totally mixed Nash equilibria of the game are those vectors 
(X\, . . . , X r ) with Xi := (xio, . . . , Xi n J for every 1 < i < r satisfying: 

(a) Xij > for i = 1, . . . , r and j = 0, . . . , nj, 

(b) Xjj = 1 for % = 1, . . . , r, 

0<j<rii 



■ ■ Xi-lj i _ 1 Xi^-\j i+1 . . . X r j r 

. . . , X r ) 
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J-i Xl 3l ■ • • x i-lji-i x i+lji+i • • • x rj r 



for i = 1, . . . , r and k = l,...,m, where 



J-i 



the sum runs over all J-i := ji ■ ■ ■ ji-iji+i ■ ■ ■ j r with < jt < n t for every t ^ i and 
(ik) ._ (i) _ W 

a j_i • C ji...j t _ 1 kj i+1 ...jr C jl-ji-lOj i+1 ...j r ' 

Let us observe that (c) is a system of n := n\ + • • • + n r multihomogeneous polyno- 
mial equations in the r groups of variables X\ , . . . , X r with m + 1, . . . , n r + 1 variables 
respectively (with degrees 1 or with respect to each group) and, therefore, it defines a 
(possibly empty) projective variety in P ni (C) x • • • x P nr (C). The complex solutions to 
the polynomial equation system (c) will be called quasi- equilibria of the game (see [7]), 
and those solutions not lying in any of the infinite hyperplanes {xio = 0} (1 < i < r) will 
be called affine quasi- equilibria of the game. 

Every quasi-equilibrium ^ := - -- ,£r) £ P ni (C) x • • • x P" r (C) determines at most 
one totally mixed Nash equilibrium of the game: for every 1 < i < r, let s& := So<j<n 4 &j 
be the sum of the coordinates of £j. If 7^ for every 1 < i < r, the unique associated 
representation of £ whose coordinates satisfy condition (b) is (Ci/sgi, • • • )£r/s§ r ); an d it 
will be a totally mixed Nash equilibrium if and only if all its coordinates are positive real 
numbers. 

2.4 On the number of solutions to a multihomogeneous system 

Let r £ N be a positive integer. Fix positive integers n\, . . . ,n r and consider r groups 
of variables Xj := (xjo, . . . , xj nj ), j = l,...,r. We say that the polynomial F G 
K[X±, . . . , X r ] is multihomogeneous of multi-degree v := (v%, . . . , v r ), where v\, . . . , v r are 
non- negative integers, if F is homogeneous of degree Vj in the group of variables Xj for 
every 1 < j < r. 

Set n := Yl T j=i n j- The classical Multihomogeneous Bezout Theorem, which follows 
from the intersection theory for divisors (see for instance |27l Chapter 4]), states that 
the set of common zeroes (over an algebraically closed field) in the projective variety 
P™ 1 x • • • x P" r of n generic multihomogeneous polynomials F%, . . . , F n with multi-degrees 
Vi := (i/ji, . . . , Ui r ) for i = 1, . . . , n is a zero-dimensional variety with 



points, where Z '■= {{ji, ■ ■ ■ ,jn) / #{k ■ jk = i} = n>i VI < i < r}. For an alternative 
proof of this result using deformation techniques, we refer the reader to [21 j . |19j and |10j . 
Note that this can be seen as a particular case of Bernstein's theorem on the number of 
common roots of sparse systems [3J . 

The quantity Bez ni .„ )TlT .(fi; . . . ; u n ) is called the Bezout number of the generic multi- 
homogeneous polynomial system. If k\,...,kt are positive integers with Y2i=i = n, we 
will use the notation Bez ni) ... )rv (z/i, k%; . . . ; Vf, kf) for the Bezout number of a multihomo- 
geneous system with ki polynomials of multi-degree Vi for every 1 < i < t. 

The equations arising in our particular setting for the computation of totally mixed 
Nash equilibria of a game are multilinear (see Section |2.3|) . Moreover, for a game with 
r players with n\ + 1, . . . , n r + 1 pure strategies respectively, we have a system of n = 



n 




(2) 
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^2j=i Tij polynomial equations consisting of exactly n, polynomials of multi-degrees equal 
to di := (1, . . . , 1, 0, 1, . . . , 1) G (No) r (where the lies in the ith coordinate of di) for every 
1 < i < r. Then, according to the previous formula, a system coming from a generic game 
with the given structure will have 

r rii 

6 --= e (nn^) (3) 

solutions in P" 1 x • ■ • x F nr , where 3 := {(in, • • • ,iim,- •• ,jrl, • • • ,jm r ) I #{jhk ■ jhk = 
i} = rii VI < % < r}. In fact, for a "generic" game, this will be the number of totally 
mixed Nash equilibria (see |18j). 

We are going to deal with the case in which 5 > 0. This inequality can be determined 
by considering the set of exponents appearing with non-zero coefficients in each of the 
polynomials in the system (see [23, Chapter IV, Proposition 2.3]) and in our particular 
case, it is equivalent to the fact that rij < J2i<k<r k^j n k = n — n j for every 1 < j < r. 
From now on, we will assume that these inequalities hold. 

Taking into account that da = for every 1 < i < r, the only n-tuples contributing 
to the sum ([3]) are those where jn~ ^ i for every 1 < k < rtj, and for each of them, the 
corresponding term of the sum equals 1. Therefore, 5 equals the cardinality of the set 

3o = {(in, • • • ,irn r ) / jik + i V 1 < k < rii and #{j hk / jhk = i} = rii V 1 < i < r}. (4) 
2.5 Geometric resolutions 

A way of representing zero-dimensional affine varieties which is widely used in computer 
algebra nowadays is a geometric resolution. This notion was first introduced in the works 
of Kronecker and Konig in the last years of the XlXth century [14J. Roughly speaking, 
it consists in a parametric description of the variety in which the parameter values range 
over the set of roots of a univariate polynomial. Now, we give the precise definition we 
are going to use. 

Let V = {£ (1) ,...,£ (<5) } CA" be a zero-dimensional variety defined by polynomials in 
x\, . . . , x n ]. Given a separating linear form £ — u\Xi + • • • + u n x n G K [x±, . . . , x n ] for 
V (that is, a linear form £ such that £(£,^) ^ £(£,^) if i ^ k), the following polynomials 
completely characterize the variety V: 

• the minimal polynomial p := ni<i<5(^ — )) £ -^Pl °f ^ over variety V 
(where T is a new variable), 

• polynomials wi, . . . , w n € K[T] with degu>j < 5 for every 1 < j < n satisfying 

V={(^( v ),...,^( v ))eK n / v eK, p( v )=0}. 

The family of univariate polynomials p, wi, ■ . ■ , w n £ K[T] is called the geometric resolu- 
tion of V (associated with the linear form £). 

In our particular setting of totally mixed Nash equilibria computation, we will not 
only deal with zero-dimensional varieties in an affine space, but we will also consider 
zero-dimensional subvarieties of P™ 1 x • • • x P""" . 
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Write £ := . . . , £ r ) G P™ 1 x • • • x with & := (&o : ■ ■ ■ : &nj (1 < * < r) to 
denote a point in the projective variety P ni x • • • x P" r . Assume that V C P™ 1 x • • • x P" r 
is a zero-dimensional variety defined by multihomogeneous polynomials in K[X±, . . . ,X r ) 
such that £jo 7^ (1 < i < r) holds for every point £ € V. Then, we may associate with V 
the following zero-dimensional variety in A n , where n := n\ + • • • + n r : 

{(Ci, • • • , C) eAV^ = (W&o, • • • , W&o) Vl<J<r,(Gn 

A geometric resolution p, W\\, . . . , w\ ni , w r \, . . . , w rrir € K[T] of this zero-dimensional 
variety will also be called a geometric resolution of V C P ni x • • • x P rar . In this case, the 
geometric resolution of V provides the following description of the variety: 

V = {[(p'(rj) : wn(rj) w lni {r])), . . . , (p'(rf) : w rl {rf) : • • • : w rrir (??)))/ rj £ K, p(rj) = 0}. 



3 On the totally mixed Nash equilibria of a generic game 

This section is devoted to the study of totally mixed Nash equilibria of generic games. In 
order to do this, we will begin by treating the payoff values as parameters and computing 
a geometric resolution of the set of quasi-equilibria of the associated generic game. This 
geometric resolution will provide a rational formula of the parameters which, for generic 
values of them, represents the totally mixed Nash equilibria of the game. 

Then, we will obtain a finite number of generic conditions ensuring that a specific 
game has the maximum possible number of totally mixed Nash equilibria. Under these 
conditions, the Nash equilibria of the game can be described by the geometric resolution 
obtained by substituting the given payoffs for the parameters in our generic formulas. 



3.1 The set of quasi-equilibria of a generic game 

Here we are going to compute a geometric resolution of the set of quasi-equilibria of a 
generic game with r players with n\ + 1, . . . , n r + 1 pure strategies respectively, where 
rii € N for every 1 < i < r. 

For 1 < i < r, 1 < k < m, let := (^.^j._ l3 -. +1 ... jV )o<it<nt be a set of new 

indeterminates and 

F k ^ := xijj . . . Xi-ij i _ 1 Xi + ij i+1 . . . x T j r , (5) 

J-i 

where the sum runs over all J-i := ji ■ ■ ■ ji-iji+i ■ ■ ■ j r with < jt < i%t for every 

(i) 

t 7& i; that is, is a generic multihomogeneous polynomial of multi-degree a, := 
(1, . . . , 1, 0, 1, . . . , 1) € (No) r (where is in the i-th coordinate). 

Note that, if {C^ j r )i<i<r,o<jt<n t are new variables that are regarded as parameters 
representing the generic payoffs of the game, from the polynomials introduced in ([5]), we 
can obtain the polynomials defining the quasi-equilibria of the generic game by substituting 

Aik) _ 

J-i ■ jl—ji-lkji+i—jr jl—ji-lOji+l—jr ^ ' 
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for every 1 < i < r, 1 < k < rii and each J_j = j± ■ ■ ■ ji-iji+i ■ ■ ■ j r with < jt < fit for 
every t ^ i. Thus, we are going to compute a geometric resolution of the solution set of 
the generic polynomial system = 0, 1 < i < r, 1 < k < and then, we will make the 
substitution Q in the result to get the desired geometric resolution of the quasi-equilibria 
of the generic game. 

Our main result is the following: 

Theorem 1 Let notations be as in Section^ There is an algorithm which computes a ge- 
ometric resolution of the set of quasi-equilibria of a generic game withr players having ni+ 
1, . . . , n r +l pure strategies respectively, within complexity 0{D 2 {D+n\ . . .n r 5 log(-D)r 2 n 4 (n 3 + 
rN))), where 

D := Y n i Bez nij .„ >r ^(d ,no]di,ni)...]di,ni-l-,...;d r ,n r ) (n := 1, d := (1, . . . , 1)), 

0<i<r 

5 := Bez ni) ... )rv (di,m;...;d r ,n r ), 
n := Y n ii 

l<i<r 

N : = n+l+ Y n i (ni + l)...(n i _i + l)(n m + l)...(n r + l). 

l<j<r 

The algorithm obtains polynomials P(T), Wij(T) G Q[C^_ V ][T] with deg T P = 5, deg T Wy < 

5 and degrees bounded by D in the parameters C-^ which are encoded by straight-line 
programs of length 0(D 2 (D + n\ . . .n r 5log(D)r 2 n (n 3 + rN))). 

(i) 

Proof. The geometric resolution of the system Ft' = 0, 1 < i < r, 1 < k < ni, will be 
obtained by means of multihomogeneous resultant computations. More precisely, we are 
going to obtain the geometric resolution associated to a linear form t which separates 
the points of the generic system by computing the resultant of our system and a generic 
linear form. The minimal polynomial of I and the polynomials Wij will be obtained by 
specializing the resultant and some of its partial derivatives in the coefficients of £. 

We introduce a set of new indeterminates := (A^, A^' : 1 < i < r, 1 < j < rij) 

which stand for the coefficients of a generic affine linear form A^ + Y i<^<r A$ Xij in 

the n = ni + • • • + n r variables x^ (1 < i < r, 1 < j < rij), and we consider the multilinear 
polynomial 

F := A^x 10 . . . x r0 + ^2 A if x w ■ ■ -Xi-ioXijXi+io . . .x r0 , 

1 < i < r 
1 < j < rii 

which is obtained by homogenizing the generic affine linear form with respect to each 
group of variables X\ , . . . , X r . 

The first step of the algorithm consists in the computation of the polynomial resultant 
11 = Res(F , . . . , F$, . . . , F[ r \ . . . , F$). According to [201 Theorem 1], this resul- 
tant coincides with the specialization of the resultant of a system of multihomogeneous 
polynomials of respective multi-degrees do = (1, ...,l),dj = (0,1,...,!) for 1 < j < 
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ni, . . . , dj = (1, . . . , 1, 0) for 1 < j < n r , in which all the coefficients of the polynomial 
of multi-degree cfo corresponding to monomials not appearing in Fq are substituted for 0. 
Thus, in order to obtain the polynomial TZ we apply an adapted version of the algorithm 
underlying the proof of |X3|, Theorem 5] (see Subsection 15. ip . With our previous notation, 
the complexity of this algorithm is of order 0(D 2 (D + n\ . . ■ n r 5log(D)r 2 n (n 3 + rN))) 
and it computes a straight-line program for TZ whose length is of the same order. 

Then, a standard procedure enables us to compute a geometric resolution of the zero- 
dimensional variety defined by the system = (1 < i < r, 1 < k < n,i) from the 
polynomial TZ: the polynomials giving the parametrization of the points in the variety are 
obtained from the partial derivatives 

, ^ dTZ , 

n ° := and Ui3 := ( - * - r ' - 3 ~ ni) ' 

dA\ dA\/ 

which we compute from the straight-line program representing TZ within the same com- 
plexity order as in the first step. These partial derivatives are encoded by straight-line 
programs of length of the same order as for the straight-line program encoding TZ (see [5]). 

Note that deg .(o) TZij < 5 for every 1 < i < r, 1 < j < n{. 
A o 

Let L :=^2 i<i< r L^Xij be a generic linear form in the variables Xij, where Lij are new 

l<j<n i 

variables, and let T be another new variable. Consider the polynomial Pl € Q[A^ ik \Lij] [T] 
obtained by specializing 

4 0) -> T, A§> ^ -L^ (1 < i < r, 1 < j < m) (7) 

in TZ. Since the resultant TZ is in the ideal (F ,F^ ] : 1 < k < r, 1 < i < n k ), then 
substituting L for T in Pl, we obtain a polynomial V which is in the ideal (F fc W : 1 < k < 
r,l < i < rik). Taking into account that deg- r (P/ / ) = 5, which is the number of solutions 
of the system = 0, it follows that Pl is a multiple by a nonzero factor in Q[A^} of 
the minimal polynomial of L. On the other hand, taking the partial derivative of V with 
respect to the variable L^, we get 

W~ = (L ' - Li *> ^ + 5) (L ' ~ L ^ Aiik) ^ e ^ : 1 < ^ < r, 1 < i < n k ). 

We conclude that making the substitution ([7|) in TZo and TZij (1 < i < r, 1 < j < m) 
we obtain polynomials which complete the geometric resolution of the variety defined by 
F^' = with respect to the generic linear form L -=Y1 i<><'- LijXij. 

Finally, in order to obtain a geometric resolution for the zero set of the generic system, 
we choose a separating linear form and we substitute its coefficients for the parameters . 
As the multihomogeneous system F^' = is generic, it has no zeroes in the hyperplanes 
xm = (1 < i < r) and we can consider its zeroes as affine points by setting Xio = 1 
(1 < i < r). Let us show that the linear form I := J2i<i< r x n separates these affine 
points. To this end, it is enough to show the existence of a choice of coefficients such 
that the induced polynomial system has the maximum number of affine roots and that I 
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(i) 

separates those roots. To see this, choose coefficient vectors for the polynomials so 
as to obtain a specific system fjf' with the maximum number of affine solutions, and take 
a linear form I € Q[xij] 1 < i < r, 1 < j < n*] separating the common affine roots of the 
polynomials \ Now, making a linear change of variables in each group (1 < % < r), 
the linear form / maps to I and the specific system considered leads to a polynomial system 
of the same structure in the new variables, which is the particular system we were looking 
for. 

Hence, the algorithm proceeds to specialize 

Af> i — ^ T, A$ h+ -1 (1 < * < r), Af) i-» (1 < i < r, 2 < j < m) 
in the polynomials 7Z, TZq and TZij (1 < i < r, 1 < j < m) to obtain new polynomials 

(i) 

giving the geometric resolution of the set of common zeros of the polynomials Ft in 
P ni x • • • x P n '\ 

If (C-^ j r )i<i<r,o<j t <n t are the parameters representing the generic payoffs of a game 
with the given structure, in order to obtain the geometric resolution of the set of its quasi- 
equilibria it suffices to substitute Cy ■ ■ —C6 a n» ~ for the coefficient -A, 

1 3l—3i-l K Ji+l—Jr J1---J1— lUji+l...Jr J-i 

for every 1 < i < r, 1 < k < n,i and each J_j = j\ ■ ■ -ji-iji+i ■ ■ ■ j r with < jt < n% for 
every t^i. In this way, the algorithm obtains polynomials P(T), dP/dT(T) and Wij(T) 
in Q[cj* j r ]\T] such that the set of quasi-equilibria of the generic game in P ni x • • • x P™ r 
is represented as follows: 

dP dP — 

{ ((^(t) : Wn(t) ■ ■ : W lni (t)), . . . , (— (t) : W rl {t) : • • • : W r „ r (t))) / i € K, P(t) = 0}, 

(8) 

where K := Q(cj*' 1 ^). The polynomials P, dP/dT and Wjj are encoded by straight- 
line programs of length 0(D 2 (D +n± . . . n r 51og(D)r 2 n 4 (ra 3 + rN))) over Q, which is also 
the order of complexity of the whole computation. The upper bounds deg r P < 5 and 
deg T Wij < 5 follow from the stated upper bounds for the degrees of 7Z and IZij in the 

variable ^4q ^. □ 

3.2 Games with the maximum number of totally mixed Nash equilibria 

The existence of games with the maximum possible number of totally mixed Nash equilib- 
ria, namely the multihomogeneous Bezout number 5 of the associated polynomial equation 
system, was proved in [18]. However, no characterization of those games has been provided. 
In this subsection, we will obtain a finite family of polynomial conditions (inequalities over 
the reals) ensuring that a given game satisfying those conditions has 5 totally mixed Nash 
equilibria. We keep our previous assumptions and notation. 

Let P(T) and Wij(T) (1 < i < r, 1 < j < Hi) in Q[cj* } -J [T] be polynomials, as in The- 
orem [U which give a geometric resolution of the set of quasi-equilibria (see ([5D) of a game 
with r players with n\ + 1, . . . , n r + 1 strategies and generic payoffs (C^ j r )i<i<r,o<j t <n t - 

Consider a specific choice of payoff values c := (c^* j r )i<i<r, o<j t <n t over R and assume 
that the polynomials P(c)(T) and Wij{c)(T) obtained from P(T) and Wy(T) by spe- 
cializing the parameters at c provide a geometric resolution of the set of quasi-equilibria 
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of the game with the given payoffs. Then, the arguments in the last paragraph of sub- 
section 12.31 imply that the totally mixed Nash equilibria of the game are those points 
(£1, . . . , £ r ) G R ni+1 x • • • x R nr+1 of the form 
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( P'(c)(t) Wg(c)(t) W ini (c)(*h n< . <r1 
U(c)(t)' #(c)(t) Si(c)(i) / 



where P' := dP/dT and Si := P' + ^i<j< ni WJy, having all their coordinates real and 
positive; that is, with t belonging to 

{t G R : P(c)(t) = 0, P'(c)(t) > 0, Wy(c)(t) > V 1 < * < r, 1 < j < m) 

or 

{iGK: P(c)(t) = 0, P'(c)O) < 0, Wij(c){t) < V '1 < % < r,l < j < rii}. 
Equivalently, t must belong to the intersection 

p| {t£R:P(c)(t)=0,(P'(c)W ij (c))(t)>0}. (9) 

l<i<r, l<j<rii 

Thus, our polynomial conditions on the payoff vector c which imply that the associated 
game has 5 totally mixed Nash equilibria will state, on the one hand, that the generic 
geometric resolution can be specialized into c leading to a geometric resolution of the set 
of quasi-equilibria of the game with the given payoffs and, on the other hand, that the 
cardinality of the set introduced in Q equals 5. 

Note that each of the sets appearing in the intersection defined in ([9]) is described by 
means of one equality and one inequality of univariate polynomials over R. In order to 
estimate the cardinality of one of these sets, we will apply the following well-known result 
due to Hermite (see, for example, [1, Theorem 4.13]): 

Proposition 2 Let p,q G R[T] be polynomials and assume that p is square-free. Let S q : 
R[T]/0>) x ^[P]/(p) — >■ R &e the symmetric bilinear map defined by S q (f, g) = Trace(My gs ), 
where for every H G R[P], Mjj ■ M.[T]/(p) — > R[T]/(p) denotes the linear map Mfj[f) = 
H-f. Then 

Signature^) = G R / p(rj) = 0, q(n) > 0} - #{r? G R / p(n) = 0, q(n) < 0}. 

The main result of this section is the following: 

Theorem 3 Under the previous assumptions and notations, there is a family of n5 + 1 
polynomials So, s\f , 1 < i < r, 1 < j < n { and \<h<5, in Q[Cj^ .j r ]i<i<r,o<jt<n t \ {°} 

with total degrees bounded by A5 2 D, such that for every vector c := (c® j r )i<i<r,o<jt<n t 
with real coordinates satisfying the conditions 

S (c) + 0, S^\c) > (1 < i < r, 1 < j < m, 1 < h < S), 

the game with r players with n\ + 1, . . . , n r + 1 pure strategies and payoff values given by 
c has 5 totally mixed Nash equilibria. 
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The polynomials Sq and can be obtained algorithmically within complexity 0(5 2 (n5 2 + 
L)) from a straight-line program of length L encoding polynomials P, Wij as in Theorem 
CO The algorithm computes straight-line programs of length 0(5 2 (5 2 + L)) which encode 
these polynomials. 

F or every 1 < i < r, 1 < k < n h let Gf G Q[CjJ_ . jr ][Xi, . . . , X r ] be the polynomial 

obtained after Fjp (see ([5])) by means of the substitution stated in ©. The system 

Gf = (1 < i < r, 1 < k < m) defines the quasi-equilibria of a generic game with r 
players and m + 1, . . . , n r + 1 pure strategies respectively. 

Let P(T) and Wy (T) (1 < i < r, 1 < j < Hi) in Q[C^.. jr ][T] be polynomials as in 

Theorem [1] which give a geometric resolution of the set of common zeros of G k over an 

algebraic closure of K := Q(C-^ 

Denote P' the derivative of the polynomial P with respect to its main variable T. Let 
us consider the resultant 

S := Re Ss ,s-i(P(T),P\T)) G Q[cg.. ip ] 

of P(P) and P'(T) regarded as polynomials in the single variable T with coefficients 
in Q[C® -J, where the subindices indicate the degrees in the variable T of P and P' 
respectively. 

Let us observe that, for every real vector c = (c^ J ^)i<^< ri o<jt <n t with S*o(c) 7^ 0, the 
polynomial P(c)(T) obtained by specializing all the coefficients of P{T) at c is a non-zero 
square-free polynomial of degree 5. Furthermore, the solution set of the system G? 1 (c) = 0, 
1 < i < r, 1 < k < n{, is a zero-dimensional sub- variety of P" 1 x • • • x P" r with 5 distinct 
points and the polynomials P(c)(T), Wij(c)(T) (1 < i < r, 1 < j < rij) give a geometric 
resolution of this variety. 

Therefore, if the condition So(c) 7^ holds, taking into account that deg T (P(c)(T)) = 
5, the arguments at the beginning of this subsection imply that the game with payoff 
vector c will have 5 distinct totally mixed Nash equilibria if and only if every root t of 
P(c) is real and satisfies (P'(c)Wij(c))(t) > for every 1 < i < r, 1 < j < rii (see ©). 
This additional condition can be restated as 

#{t G R : P(c)(t) = 0, (P'Wij(c))(t) > 0} = 5 V 1 < i < r, 1 < j < m. 

Since deg T P = 5, this last condition is equivalent to 

Signature^p/v^c)) = 5 V 1 < i < r, 1 < j < ni. (10) 

Note that for any fixed pair of indices being Spiy^..^ a bilinear form defined over 
a vector space of dimension 8, the equality Signature^p/iy.^)) = 5 holds if and only if all 
the eigenvalues of the matrix of Spiyy..^ in an arbitrary basis of the space are positive, 
which is in turn equivalent to the fact that the coefficient sequence of the characteristic 
polynomial of this matrix, whose roots are all real, has 5 sign changes. The conditions 
stating these sign changes will be the inequalities ensuring that the game has the maximum 
number of totally mixed Nash equilibria. 
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Now, we will obtain the polynomials sff giving the inequalities in the statement of the 
theorem by considering the bilinear forms (Sp'w^ f° r generic payoffs and performing all our 

computations over the parameter field K = Q(cj^ that is, we will consider the bilinear 
form S P > Wi . : K[T]/(P) x K[T]/(P) -► K. Let £ K 5xS be the matrix of S P , Wij in 

the basis B := {1,T, . . . , T 5-1 }. Note that, for every real vector c = (c^ j r )i<«<r, o<j t <n t 
with S*o(c) 7^ 0, the coefficient of T" 5 of the polynomial P(c)(T) is not zero. Then, the 
matrix of the bilinear form «Sp/^ .( c ) can be computed by specializing M.^ at c. 

Actually, to deal with polynomials instead of rational functions of the parameters, we 
will not compute the matrices of the bilinear forms Spiy/ i . but (scalar) multiples 

of them with polynomial entries. In what follows, we will first show how to obtain these 
matrices. Then, we will compute the coefficients of their characteristic polynomials, and 
we will estimate their degrees and the overall complexity for their computation. 

Before proceeding, we are going to make some remarks about the denominator-free 
computation of traces of multiplication maps in K[T]/(P). It is easy to see that for all 
H € K[T], we have M# = H(Mt)\ that is, every multiplication map is a polynomial 
function of the multiplication map Mr : K[T]/ (P) — > K[T]/(P), whose matrix in the basis 
B := {1, T, . . . , T 6 - 1 } can be read off from P. In fact, if we write P = ELo^fc T * : > with 
p k € Q[Cj l ^ jr ]i<i< r ,o<j t <nt, tlie matrix of M T in the basis B is 



/ 
1 





V o 



-2a \ 



PS 

EL 

PS 




1 



PS-l 

ps 



J 



In order not to deal with denominators in our computations in the case when H G 







][T], we will consider the multiplication map M Ps t, whose matrix in the basis 



B is 



M := 



( o 

PS 





-Po 
-Pi 



\ 





Ps 



-Ps-i J 



VH = YLo h hT h with 6 fc €Q[C (i) 



ji...>j 



we have that Trace(M#) = Y? h=() b h Ttace({M T ) h ) 
and, therefore, we can obtain a multiple of Trace(M#) working over the polynomial ring 

Q[cj^ j r ] due to the following identity which holds for every K > d: 



pfTrace(M H ) 



d 



K- 



h=0 



h b h Ti a ce((M nT ) h ). 



(11) 



Taking into account the previous remark, we are able to obtain a multiple of the matrix 
by a sufficiently big power of the leading coefficient ps with all its entries in Q[cj^ -J 
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by working over this polynomial ring: for every 1 < a, (3 < 5, we have that 

(M^) a ,p = Trace(M TQ _ lp , H/8jT ,- 1 ) (12) 

and then, as deg T (P') < 6 and deg T (Wij) < 5, the upper bound deg(T a - 1 P , W i jTP- 1 ) < 
45 — 4 holds. We conclude that the matrix p^~ 4 A4^ 1 ^ has all its entries in the polynomial 

ring Q[Cj^ -J. Moreover, as we are multiplying the matrix JVi^fi by an even power of pg, 
for every payoff vector c with S'q(c) ^ 0, the identity 

Signature(.M (y) (c)) = Signature^c) 4 ' 5 - 4 .^^)) 

holds (note that Sq(c) / implies pg(c) / 0). 

Therefore, our algorithm will compute the matrices pf~ A and then, their char- 
acteristic polynomials X i5 = T 5 + Y?h=i(- l ) h S \f T& ~ h € ®[ C fl..j r \ \ T \ for ever y 1 < * < 
r, 1 < j < rij. Our previous arguments imply that for a payoff vector c with Sq(c) ^ 0, 
condition (fTUj) is equivalent to 

^(c) > V 1 < i < r, 1 < j < m, 1 < h < 5. 

Now we detail the successive steps of the algorithm, we estimate its complexity and 
the length of a straight-line program encoding the output polynomials sf- \ and we give 
an upper bound for the degrees of these polynomials. 

The entries of the matrices p^~ 4 A4^ are obtained according to identities (|12l) and 

•CD. 

First, from a straight-line program of length L encoding the polynomials P and Wij 
(1 < i < r, 1 < j < rtj), we obtain a straight-line program of length 0(5 2 L) for the 
coefficients of P in the variable T and a straight-line program of length 0(5 2 (L + n)) for 
the coefficients of all the polynomials P'Wij (1 < i < r, 1 < j < n,) in the variable T 
within complexity 0(5 2 (L + n)) (see [U Lema 21.25]). 

Then, using the coefficients of P, the matrix J\A 6 Q[cj^ jj' 5 *' 5 °^ ^ e mur tiplication 
map M Ps t is constructed and the powers M h for h = 0, . . . , 45 - 4 are computed. Note 
that every matrix in Q[c!j^ jr ] 5x5 can be multiplied by Ai with 0(5 2 ) additions and 

multiplications in Q[cj^ -J and so, the powers .M' 1 for /i < 45 — 4 of this matrix can 
be computed within 0(5 3 ) polynomial additions and multiplications. Thus, we obtain 
straight-line programs of length 0(5 2 (5 + L)) for the entries of the matrices M (0 < h < 
45 — 4) within complexity 0(5 3 ). 

Therefore, a straight-line program of length 0(5 2 (5 + L)) computing Trace ((M ps T) h ), 
for every < h < 45 — 4 is obtained with 0(5 2 ) additional operations. 

Now, fix i,j with 1 < i < r and 1 < j < m. Let P'Wij = Y?h=o $ S) 'T h be the 
expansion of P'Wij in the variable T with coefficients in Q[cj^ ^J. Recall that in the first 

step of the algorithm we have obtained a straight-line program of length 0(5 2 (L + n)) for 

(a) 

the polynomials oi . For every a, (5 with 1 < a, (3 < 5, we obtain a straight-line program 
for (p 45_4 A / I^^)a, / 3 by means of the formula: 

2<5+a+/3-4 

pf- 4 Trace(M r «- lpWi . r ,- 1 ) = £ pf- 4 " ?t 6£ ) a _ /3+2 Trace(M feT)h ) 

h=a+l3-2 
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within additional complexity 0(5). The length of this straight-line program is 0(5 2 (5 + 
L + n)). 

The complexity of the computation of the entries of all the matrices p$ 5 ~ 4 \A4^ from 

the coefficients b£ and the traces Trace((M Pi5 ;r) /l ) is 0(nd 2 ). 

Finally, we apply the division-free algorithm described in [2] in order to compute the 
coefficients of the characteristic polynomial Xy of each of these matrices. For every 1 < 
i < r, 1 < j < rii, the algorithm produces straight-line programs of length 0(5 2 (5 2 +L+n)) 
encoding these coefficients within complexity 0(<5 4 ). 

The overall complexity of the algorithm is 0(5 2 (n5 2 + L)). 

Let us observe that the entries of the matrices p^ s ~ i A4^^ are polynomials in Q[cj^ -J 

of total degrees bounded by 45D and so, the polynomials S^j obtained by the described 
procedure, which are the (signed) coefficients of the characteristic polynomials Xy, have 
total degrees bounded by A5 2 D. □ 

Note that, with the same notation as in the previous theorem, for a generic game 
with the given structure (namely, any game whose payoff vector c satisfies So(c) ^ 0) 
the conditions S^\c) > are equivalent to the fact that the considered game has the 
maximum possible number of totally mixed Nash equilibria. 

4 The set of totally mixed Nash equilibria of an arbitrary 
game 

When dealing with a particular game with specific payoff values, the geometric resolution of 
the set of quasi-equilibria of a generic game with the same structure might not be useful 
in order to describe the quasi-equilibria of the given game. For instance, the minimal 
polynomial P(T) or its derivative ^(T) computed by the algorithm shown in Theorem [TJ 
could vanish identically when specialized at the considered payoff values. This section is 
aimed at adapting the procedures previously developed in order to handle these particular 
situations. 

First, we show an algorithm computing a geometric resolution of the set of affine quasi- 
equilibria of the game provided that it has a zero-dimensional set of quasi-equilibria. In 
a second step, this algorithm is extended to find all isolated affine quasi-equilibria of an 
arbitrary game. 

Once a finite set of affine quasi-equilibria of a game including the isolated ones is 
computed, we apply standard algorithms of semialgebraic geometry to estimate the number 
of isolated totally mixed Nash equilibria of the game. 

4.1 Games with a zero-dimensional set of quasi-equilibria 

As in the previous sections, we consider an r-person non-cooperative game in normal form 
in which the players have n\ + 1, . . . ,n r + 1 distinct available pure strategies each. Let 
cM := [cfV j r )o<j k <n k denote the payoff matrix to player i for every 1 < i < r. 

Then, the polynomial equations defining the Nash equilibria of the game (see Section 
I2.3P can be obtained by specializing the coefficients of each of the generic multilinear 
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polynomials F£' introduced in ([5]) in the vector := ( a f^j i _ 1 j il j r )o<j t <n t defined 



(ik) ._ Ji) Ji) 



Jl.. .ji-! kji+1 . . .j r + 1 ■ ■ -j 



from the payoffs of the game as follows a - ■ ■ ■ := c 

Thus, if a := (a^)i<i< r , i<k<m, the set of quasi-equilibria of the game is the algebraic 
variety 

V a := {£ := • • • , £r) € P™ 1 x • • • x P""- / F fc (i) (a,e) = 0Vl<i<r, l<fc< nj. 

In this section we consider the case when V a is zero-dimensional. In order to decide 
whether this is the case for given payoff values, we study the non vanishing of an adequate 
multihomogeneous resultant. More precisely, we consider a generic polynomial of multi- 
degree do := (1, . . . , 1) in the groups of variables Xi, . . . , X r , 

Fo= Yl A fi.jr x ih ■ ■ ■ x rj r 

1 < i < r 
< ji < n t 

and the resultant R := Res(F , F^\. . . , , . . . , F x (r) , . . . , Let R a (Af^ jr ) be the 

polynomial obtained by substituting the coordinates of a = (a^ fc ^)i<j< r , i<fc<n 4 for the 
variables in the polynomial R. Then: 

Remark 4 The algebraic variety V a is either zero- dimensional or empty if and only if the 
polynomial R a is not identically zero. 

Proof. If V a is empty, the result is straightforward. If V a is zero-dimensional, there exists 
a multilinear polynomial /q € Q[-Xj., . . . , X r ] which does not vanish at any of the (finitely 
many) points of V a and therefore, R a does not vanish at the coefficients of /o- 

On the other hand, if the projective variety V a has positive dimension, any multilinear 
polynomial /o has zeros in V a . Therefore, R a is identically zero. □ 

Now we will show how to decide whether R a is identically zero algorithmically and we 
will estimate the complexity of this procedure. 

First, we compute the multihomogeneous resultant R by means of the algorithm de- 
scribed in [13] adapted according to Subsection 15.11 below. Thus, a straight-line program 
of length C := 0(D 2 (D + n\ . . . n r 51og(D)r 2 n 4 (n 3 + rN))) encoding R is obtained within 
complexity of the same order as C The specialization A^ := a^' which leads us to the 
polynomial R a does not modify this complexity order. 

Since R a is a multivariate polynomial in ni<i<r( n * + l) variables encoded by a straight- 
line program, we do not know its coefficients; moreover, the complexity to compute them 
is exponential. Therefore, in order to decide whether it is the zero polynomial or not, we 
will make a suitable specialization so that solving this problem amounts to solving the 
same problem in the univariate setting. 

More precisely, we specialize the variable coefficients A^ into successive powers of 
a new variable t, that is, 

A (0) • = ^l + ("'l+ 1 )j2 + (™l+ 1 )("2+ 1 )^ + --- + (n^o 1 (™J+ 1 ))>. (13) 

3x***3r ' V / 
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Making this specialization in Fq, we obtain a new polynomial /q G . . . ,X r ]. Let 

us show that if V a consists of finitely many points, there is a polynomial /o(£o) with 
coefficients of this type (namely, successive powers of the same scalar to) which does not 
vanish at any point of V a . 

For every £ := . . . , £ r ) G P™ 1 x • • • x P" r , we have 

/o(0(t)= ^ eiii---Cr>^' 1+(ni+1)j ' 2+(ni+1)(n2+1)i3+ ''' +(n ^ eC[t], 

which is a non-zero polynomial, due to the fact that there exists at least one choice 
of ji,---,j r for which the product £1^ ■■■£ r j T is not 0. Now, if V a is a finite set, we 
may consider the polynomial A(t) := n^ev a /()(£)(*)• Note that A G C[i] is a non-zero 
polynomial and, therefore, there exists to G Q with A(to) 7^ 0. Due to the definition of A, 
the polynomial fo{to) € Q[^i, • • • ,X r ] does not vanish at any point of V a . This implies 
that the resultant of fo(to), (o^) (1 < i < r, 1 < < m) is not zero, that is, the 
polynomial R a ,t G Q[t] obtained by specializing i? a following (fT3|) does not vanish at to- 
We conclude that R tt) t is not the zero polynomial. 

Therefore, in order to decide whether R a G Q[A^ .J is zero or not, it suffices to decide 
whether the associated polynomial R a j G Q[t] is zero or not. Taking into account that 
deg(R a ,t) < (lli<j<r( n * + 1) — 1)^) ^ follows that this task can be achieved by evaluating 
a t (lli<i<r( n * + 1))^ distinct values of t, which is done by substituting the powers of 
these values for the variables ■ according to (|13h in the straight-line program for R a . 



3l—3r 

The overall complexity of this procedure is ^(rii<i<T-( n * + -0) Jr ^) (Tli<i<r( n i + 1))^ = 
0(D 2 (m . . . n r ) 2 5(D + m...n r d log(L>)r 2 n 4 (n 3 + rN))). 

Assume now that the set of quasi-equilibria of the given game is finite, that is, the 
polynomial R a is not identically zero. 

Proposition 5 Following the previous assumptions and notation, there is an algorithm 
which computes a geometric resolution of the set of affine quasi-equilibria of a game with r 
players having n\ + 1, . . . , n r + 1 pure strategies respectively within complexity 0(5 s D 2 (D-\- 
m . . . n r 5log(D)r 2 n 5 (n 3 + rN))) provided that the associated set of quasi-equilibria V a is 
zero- dimensional. 

Proof. We will obtain the desired geometric resolution from the specialized resultant R a . 
In order to do this, we first show that this specialized resultant factorizes as follows: 

R a = C a \{ F (er«, with C a G Q, m s G N. 

Set T a ■= ri£ev -^b(£) £ Q[^°^]- For a given coefficient vector a^°\ we have that 
•^a(a^) = if and only if there is a point £ G V a such that Fo(a(°),£) = 0; that is, 
if and only if F (at®), F^ 1] (a^), . . . , F£ ) (a< rnr >) have a common root in P ni x • • • x P™'\ 
But this last condition is equivalent to the fact that the resultant Res(i ? o, F^\ . . . , F^ ) 
vanishes at a(°\ a^ n \ . . . , a^ rrir ^ or, equivalently, that R a (a^) = 0. We conclude that the 
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polynomials T a and R a have the same zero set and so, they have the same irreducible 
factors, which proves our assertion. 

We are interested in describing the set of affine quasi-equilibria of the game, that is, 
the variety '■= {£ S K : («o / V 1 < i < r}. For this reason, we are going to compute 
the polynomial Rf := H^y^e F (£) m t. Note that 

R a = C a H F (£H [J F ^) m5 , (14) 

where = n§eV aff -Fb(£) m? is monic in the variable A °/ an d T[^v a -V* a ^b(£) m? does 
not depend on this variable. Then, C a Y\^ &Va _ Va ,s i<b(£) me is the leading coefficient of R a 

in the variable A °\. Therefore, R^ can be obtained from R a by dividing it by this 
leading coefficient. 

Let us consider the generic form fo := A^xio ■ ■ ■ %r0+ A^yxio . . . o^ij^i+i o • • • %r0- 



1 <i < r 

i < j < m 



After specializing as follows: 

4 (0) ,(0) 

Aq qjq q i— > where the index j is in the ith place (1 < j < nj) (15) 

• i y otherwise 

we obtain the polynomial P a := ngey aff fo(0 ms - The geometric resolution of the variety 

can be obtained from the square-free part P a := W^y^s /o(£) °f an d the partial 
derivatives of P a as shown in the proof of Theorem [1] and substituting the powers of a 
conveniently chosen scalar for the variables Aij. 

Now we estimate the complexity of the different steps of our algorithm. 

Computation of R^: First, we compute the degree d a := deg (o) (R a )- To do this, let 

A o...o 

to € Q be the element we have already obtained such that R a ,t (to) 7^ 0. Let Ra °^ € Q[v!q ^ ] 
be the (non-zero) polynomial obtained from R a after specializing it as in (|13p for every 

(ill ■■■ ijr) 7^ (0, ... ,0) setting t = to. Because of JHJ), deg(J2a ) = d a . Then, in order 
to compute d a , it suffices to compute the coefficients of R a to ^ up to degree 5 (which is 
an a priori upper bound for d a ). The complexity of this computation is of order 0{5 2 £). 
Now we obtain an sip of length 0(5 2 C) for the coefficient of (A^ ) da of R a and we obtain 
Ra S by dividing R a by this coefficient. As the divisor does not vanish when its variables 
are specialized in the successive powers of to, this polynomial division can be done within 
complexity 0(5 4 C), and produces an sip of the same order ([28]). 

Computation of the polynomial P a : We obtain P a as the quotient P a /gcd(P a , dP , a m ) . 

3 

The specialization in f)15 j) gives an sip for P a without modifying the previous length. 

To obtain the square- free part of P a , we apply a well-known subresultant-based pro- 
cedure for the computation of the gcd of two polynomials (see, for instance, [1]). We 
proceed as before by specializing all the variables in P a , except for the main variable A °^ 
into the successive powers of a new variable t and deal with the polynomial we obtain in 
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Q[t]L4[, 0) ]. We will use the following notation: given F G Q[Ajp,A$], F® will denote 
the polynomial obtained from F by doing the previously mentioned specialization. 

Let G = gcd(P a , As, for & ^ £ 2 we have that /^(fi) y fo\&) are relatively 

prime irreducible polynomials in Q[^4i , t], we have that 



8A {0) 



G(t) = II fo (er £_1 =gcd(pi 



Then, we can obtain d a := deg(gcd(P a , |fo ) ) using the polynomials P a and — fey. In 

8A dA 

order to compute the degree of the gcd of these polynomials, we look for their first non- 
zero subresultant. In each step, to decide whether the considered subresultant (which is a 
polynomial of degree at most 25 2 n in Q[t]) is zero or not, we evaluate the variable i in a 
sufficient number of elements of Q. First, we obtain an sip of length 0(5 4 C+n) for Pa" 1 and 
then an sip of length 0(5 2 (5 4 C + n)) for its coefficients in the variable A^. For a specific 
evaluation of t the complexity of the computation of all the subresultants is 0(5 e C + 5 2 n), 
and therefore, the whole complexity of this step is bounded by 0(5 8 Cn + <5 4 n 2 ) . We obtain 
G = sr a . gcd(P a , ® p £s ) as the polynomial subresultant of index d a , where sr a is the main 

subresultant of index d a , within complexity 0(<5 6 £). 

Finally, the polynomial P a is obtained by dividing sr a P a by G. Note that we already 
know a point r where sr a does not vanish (r is the value obtained in the previous step 
when computing d a ). Then, if we consider the nonzero polynomial G^(Aq,t) € Q[Aq ^] 
of degree d a , by evaluating it in at most d a + 1 elements in Q, we obtain £ Q such that 
GW(o,q , t) 7^ 0. This enables us to compute the quotient P a by applying the classical 
division avoiding algorithm in [28]. The complexity of this step is of order 0(5 S C). 

Observe that the linear form £ whose coefficients are (— T fc )o<fc< n _i is a separating 
linear form for V^ aff . To compute the required geometric resolution associated with £ from 
the polynomial P a we proceed as in the proof of Theorem [TJ we specialize P a and its 
partial derivatives at the coefficients of £. 

The overall complexity of the algorithm is 0(5 s Cn+5 4: n 2 ) = 0(5 8 D 2 (D+n\ . . . n r 5log(D)r 2 n 5 (n 3 + 
rN))). □ 

Now, we are able to compute the number of totally mixed Nash equilibria of a game 
with zero-dimensional set of quasi-equilibria from the geometric resolution given by Propo- 
sition [5j 

Let p and Wij (1 < i < r, 1 < j < rii) be the polynomials giving the geometric 
resolution of V^ s . Then, the totally mixed Nash equilibria of the game are the points 
(£i, ...,£r) with & = (p'(t)/si(t),Wii{t)/si(t), . . . ,w ir (t)/si(t)), where Si = p' + YJjL\ w ij 
for i = 1, . . . , r, and t is a root of p, having all their coordinates real and positive. 

Then, the number of totally mixed Nash equilibria of the game equals the cardinality 
of the union of the sets 

{t£R: p(t) = 0,p'(t) > 0, Wij(t) > VI < % < r, 1 < j < rii) 
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and 

{t € R : p(t) = 0,p'(t) < 0, Wij(t) <0\/l<i<r,l<j< nj. 

Once we have the polynomials p, p' and w^j (1 < i < r, 1 < j ' < rtj}) encoded in dense 
form, by using the algorithm in [6, Section 3.3], it is possible to compute this cardinality 
within complexity 0(n5 3 ). 

Therefore we have: 

Proposition 6 Following the previous assumptions and notation, there is an algorithm 
which computes the number of totally mixed Nash equilibria of a game with r players having 
ni+1, . . . , n r +l pure strategies respectively within complexity 0(5 s D 2 (D+ni . . . n r <51og(D)r 2 n 5 (n 3 + 
rN))) provided that the associated set of quasi- equilibria of the game is zero-dimensional. 

4.2 Computing the isolated affine quasi-equilibria of an arbitrary game 

When we are dealing with an arbitrary game, it may happen that the set V a of its quasi- 
equilibria has positive dimension and, therefore, the polynomial R a introduced in the 
previous subsection is identically zero. In this case, we are going to use a procedure 
applied in [TO] to obtain a non-zero multiple of the minimal polynomial of a generic linear 
form over the isolated points of the set of affine quasi-equilibria. For the sake of 
completeness, we are going to explain this procedure briefly. 

We consider a sufficiently generic coefficient vector b := (6^ fc ^)i<j< ri i<fc<n^ such that 
Rb (we remark that the vector b can be chosen at random or effectively constructed 

as the coefficient vector a system with 5 many common roots). Then, if r is a new 

(i) 

variable, the polynomial R a +r(b-a) i s non-zero. If denotes the polynomial obtained 
from by evaluating x j0 = 1 (1 < j < r), let ff ] := f£\aW + r(ft( ifc ) - a^)) 
(1 < i < r, 1 < k < rtj). Let P a + T (b-a) be the polynomial obtained from R a+T ^- a ) by 
specializing it as in (fT5|) and let L := Ylij ^ij x ij- 

As P a + T (b~a) ( r ) AW) | t =o= 0, there exists e > such that P a + T (b~a) = T P f° r a 
polynomial P € Q[r, A^] with P(r, A^) | T=0 ^ 0. Let P := P I _ n . m _ „. We are going 

to show that this polynomial P is a multiple of the minimal polynomial of L over the set 
of isolated points of . 

Since the polynomial R is a linear combination of Fq, Ffp (1 < i < r, 1 < k < nj), we 
have that P a+T{b _ a) \ A Q) = _ L = T e P \ A (o) = _ L & {Jl ■ 1 < i < r, 1 < k < th) G Q[t, Xij : 1 < 
i < r, 1 < j < rii] . 

Now, taking into account that the number of variables r,Xij (1 < i < r, 1 < j < nj) 
is greater than the number of generators of the ideal (fj?' 1 : 1 < i < r, 1 < k < n«), we 
deduce that each irreducible component of the variety V T defined by this ideal in A n+1 
has dimension at least 1. Then, for each isolated point £ of V^ aff , the point (0,£) € A n+1 
lies in an irreducible component C of V T such that r ^ P\C). Therefore P I .(o)_ € ^(C) 

and P |t=l vanishes at £. 

Now, we are going to show how to compute the polynomial P. 

Let P a + T (b-a) = Y^j=oPj{ T ){-^Q } y i where Pj(r) are polynomials in Q[^4^][t] such that 
Pj(r) ^ for some < j < D. Then e := max{/c : r k divides pj V0 < j < D}. 
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The procedure for the computation of P from a straight-line program encoding R runs 
as follows: 

(1) Obtain an sip for P a + T (b-a)- 

(2) Obtain an sip encoding the coefficients Pj(r) (0 < j < S). 

(3) For j = 0, . . . <5: 

(3. a) Obtain an sip for the coefficients of Pj(r) = Y^k=oPjk jk ■ 
(3.b) Compute ej := min{k : pjk ^ 0}. 

Compute e := minjej : < j < 5}. 

(4) Obtain an sip for P = Y^ S j=o(~ ^YPjeT^ from the sip's encoding pj e (0 < j < 8). 

Now, we are going to estimate the complexity of this procedure. Let C be the length 
of a straight-line program encoding R. 

In step (1), we compute first a + t(6 — a) with 3N operations and then, we specialize 
the sip for R according to (fT5j) . We obtain an sip of length 3./V + C for P a + T {b~a)- The 
interpolation in step (2) takes 0(5 2 (N + C)) operations. The complexity of step (3. a) is of 
order 0(D 2 5 2 (N+£)). In order to determine whether each of the multivariate polynomials 
Pjk £ encoded by an sip is zero or not, we apply the probabilistic Zippel-Schwartz 

zero test (see [M])- The overall complexity of step (3) is 0(D 2 S 3 (N + £)). Step (4) does 
not modify the order of the complexity. 

From the polynomial P we have computed we can obtain a geometric resolution for a 
finite set of points including the isolated points of V^ aff in the same way we showed in the 
proof of Theorem [lj obtaining a description for the isolated affine quasi-equilibria of the 
game. Therefore, we have: 

Theorem 7 Let notations be as in Theorem [7J There is a probabilistic algorithm which 
computes, given a game with r players having m + 1, . . . ,n r + l pure strategies respectively, 
a geometric resolution of a finite set of points including the isolated affine quasi-equilibria 
of the game within complexity 0(D 4 5 3 (D + n\ . . . n r 5log(D)r 2 n 4 (n 3 + rN))). 

Proceeding as in the previous section, from the geometric resolution given by this 
theorem, we are able to obtain an upper bound on the number of isolated (in the complex 
space) totally mixed Nash equilibria of the game within the same order of complexity as 
in the previous theorem. 

5 Resultants in the multihomogeneous setting 
5.1 Computing resultants 

In this subsection we are going to show how the algorithm in [13] can be adapted in order 
to compute the multihomogeneous resultants we use. That procedure computes multiho- 
mogeneous resultants under the assumption that the coordinates of each multidegree are 
all positive and this is not the case in our situation. 
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First, we will prove that the multihomogeneous resultant we want to compute is a 
non-constant polynomial and then, we will show how to modify the Poisson formula on 
which the algorithm in [13] relies All the computations required once the Poisson formula 
is recursively applied run in the same way as in [13j . 

In order to do this, we are going to use the theory in [30] and [20]. We can ap- 
ply these results to our situation because the multihomogeneous resultant of a family of 
multihomogeneous polynomials Go, ... , G ri in r groups of variables Xj = (xjo, . . . , xj nj ), 
with X]j=i n j = n > coincides with the sparse resultant of the dehomogenized polynomials 
go, . . . , g n obtained by setting Xjo = 1 for every 1 < j < r. 

Let Aq, ■ ■ ■ ,An C 7L n be finite sets and let go,-..,g n be polynomials with supports 
Aq, ■ ■ ■ , A n respectively. 

For any subset J C {0, . . . , n}, let Cj be the lattice generated by ^2j e jAj. Follow- 
ing [30J, if I C {0, ...,n}, the collection of supports {Ai}i^i is said to be essential if 
rank(£/) = #1 — 1 and rank(£j) > # J for each proper subset J of /. 

If there is a unique subcollection {Ai}i^i which is essential, the resultant Res(go> • • • , 9n) 
coincides with the resultant Res(gj;i € I). 

Note that if G is a multihomogeneous polynomial in r groups of n\ + 1, . . . , n r + 1 
variables each with multidegree (z^i, . . . , u r ), the set of exponent vectors of the dehomoge- 
nized polynomial g is A nijjyi x • • • x A nr Ur C Z ni x • • • x 7/ lr , where A ni}Ui is V\ times the 
unitary simplex in for i = 1, . . . , r. 

In our particular case, we will deal with families of multihomogeneous polynomials in 
r groups of m>i + 1, . . . , m r + 1 variables respectively and multidegrees (1, . . . , 1) or d% = 
(1, . . . , 0, . . . , 1) (with in the ith coordinate) for i = 1, . . . , r. If one of the polynomials 
in the family indexed by / has multidegree (1, . . . , 1), then rank(£/) = ^I=i m i> an d the 
same holds if the family contains polynomials with multidegrees di and dj with i ^ j. If 
all the polynomials in the family have multidegree di, then rank(£/) = Ylj^i m j- 

The resultant will be computed recursively by applying Poisson's formula ([20, Lemma 
13]). At each step, we will have to compute a multihomogeneous resultant in one of the 
following settings, where mi, . . . , m r G N and m := m\ + • • • + m r : 

(I) m + 1 multihomogeneous polynomials in r groups of mi + 1, . . . , m r + 1 variables: one 
polynomial with multidegree (1, . . . , 1) and, for every 1 < i < r, nii polynomials with 
multidegrees d{ := (1, . . . , 0, . . . , 1) with in the coordinate i, under the assumption 
that rrii < Ylj^i m j f° r ever y 1 ^ i < r i 

(II) m + 1 multihomogeneous polynomials with multidegrees (1, . . . , 1) in r groups of 
mi + !)•••) m r + 1 variables each, 

(III) m multihomogeneous polynomials in r groups of variables with mi,ni2 + l, ■ ■ ■ , m r + l 
variables each, with rrii polynomials of multidegree di for every 1 < i < r, under the 
assumption that rrii < Ylj^i 171 j f° r every 1 < i < r. 

Now we are going to start with the recursion. 

The resultant 1Z = Res(i ? o, F-J , . . . , Fn x , • • • , F[ J , . . . , Fn)) we want to compute in 
the proof of Theorem [1] involves a family of multihomogeneous polynomials satisfying the 
conditions in (I) with rrii '■= n i and m := n. Note that, for a generic game to have a 
non-empty set of quasi-equilibria, the inequalities rrii < ^2j^i m j f° r ever Y hold. 
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Let Gq and (1 < i < r, 1 < k < rrii) be a family of polynomials satisfying the 
conditions in (I) and let 2 := {(0, 0)} U {(i, k) : 1 < i < r, 1 < A; < m,j}. 

First, note that rank(£j) = m = j^I — 1. Let J be a proper subset of X. If there 
exist (i,k), (i',k') in J with i ^ i', then rank(£j) = m > j^J and the same holds 
if (0,0) £ J. On the other hand, if J C {(i,k) : 1 < k < rrii} for a fixed i / 0, then 
rank(£j) = Ylj^i m j — m i — #J- Therefore, the set of all supports is the unique essential 
subset. So, the resultant is not constant and the following identity holds: 

Res(GQ°' ) , (Gj^)l<i< r; l<fc< mj ) = Jj9o°' ) (£) ] [ R- es (Cij- )l<!<r,KKm,), 

l<j<r 

where <7q is the dehomogeneized polynomial obtained from by evaluating X£ me = 1 
(1 < £ < r), V is the set of common zeros in A m of the polynomials g$ obtained in the 
same way from the G^ 1 (1 < i < r,l < k < mi), and, for each 1 < j < r, is the 
polynomial obtained from by setting Xj mj = 0. 

(La) Note that if rrii > 2 for every 1 < i < r, (up to renaming variables and polynomi- 
als) each of the resultants Res(G^j) involves a family of polynomials satisfying the 
conditions in (III). 

(Lb) Without loss of generality, assume now that rri\ = 1. Here, we can discard the first 
group of variables. Then, we have to compute the resultant of m polynomials in r — 1 
groups of rri2 + 1, . . . , m r + 1 variables each with rrii polynomials with multidegree 
(1, . . . , 0, . . . , 1) where the zero is in the (i — l)th coordinate for 2 < i < r and one 
with multidegree (1, . . . , 1). 

If, for every 2 < i < r, rrii < Si<j<r,jVi m i' smce m i = 1 we deduce that rrii < 
X^2<j<r jjti m j an d so ! the polynomial system obtained is of the form (I) but with 
one group of variables less than the original one. 

On the other hand, if rrii = 1 + S2<j<r j^i m i f° r some 2 < i < r, we have that 
rrij < rrii for every j ^ i. Therefore, the unique essential subset is {(i, k) : 1 < k < 
rrii} and the resultant to be computed is the resultant of the corresponding family 
of rrii polynomials of multidegree (1, . . . , 1) in r — 2 groups of rri2 + 1, . . . , rrii-i + 
1, mj_|_i + 1, . . . , m r + 1 variables each, which is the situation in (II). 

If the polynomials involved satisfy the conditions in (II), all the coordinates of their 
multidegrees are not zero and so, we can apply the algorithm in |13j for the computation 
of their resultant. 

To analyze (III), let us consider first the case when r = 2. Here, the assumption on 
the numbers rrii implies that rri\ = m<i := M. 

(III. a) We consider the resultant of M polynomials with multidegrees (0, 1) and M polyno- 
mials with multidegrees (1,0) in two groups of M and M + 1 variables respectively. 
Now the unique essential set is the corresponding to the first M polynomials and, 
therefore their resultant equals the determinant of the matrix of their coefficients. 
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Assume now that r > 2. Note that the equality rrii = m j ma y be valid for at most 
one value i. If, on the contrary, = ^j^n mj and mj 2 = X^yi 2 m i n °ld for ii 7^ 12, it 

follows that *}2j J Li 1 i 2 mj = 0, which implies r = 2. Let (l<i<r, l<fc< mj) be a 
family of polynomials satisfying the conditions in (III). 

(Ill.b) If mi = 1, we are under the same assumptions as in (Lb). 

(III.c) If mi > 2 and mj = Ylj^i m j f° r some 2 < % < r, the set {(i, k) : 1 < k < rrii} 
is essential. The assumption r > 2 implies that it is the unique essential subset. 
Then, the resultant we compute is that of rrii multilinear polynomials of multidegrees 
(1, . . . , 1) in r — 1 groups of mi, m.2 + 1, . . . , mj_i + 1, m^+i + 1 . . . , m r + 1 variables 
respectively. Thus, we are in situation (II). 

(Ill.d) If mi > 2 and rrii < Ylj^i m j f° r every 2 < i < r, we deduce that rrii < ^1 — 1 + 
Sjyi j m i f° r every 2 < i < r. Therefore, the unique essential subset is the whole 
family of polynomials and applying Poisson's formula we obtain: 

l<i<r;l<k<nii 

)=n^i (o n Res (( G i/ )2<fe<m 1 ;(Gg ) ) 2 < i < r ,i< fe < mi ), 

^eV 2<i<r 

where 5^ is the dehomogeneized of by setting Xj mj = 1 for every 1 < j < r, 

V is the set of common zeros in A m_1 of the polynomials (2 < k < mi), 
(2<i<r,l<k< rrii) obtained in the same way from the homogeneous polynomias 

(i) (i) (i) 

G k , and G k { is the polynomial obtained from G k by setting xi„ H = 0. 

According to the previous formula, we will have to compute r — 1 multihomogeneous 
resultants of rrr! = m — 1 polynomials each: 

For I = 2, . . . , r, setting m! := m — 1, m' x := mi — 1 and m' i := rrii for i ^ 1, the 
resultant to be computed involves rrr! polynomials in r groups of m! x + 1, . . . , m! l _ x + 
1, m,;, m' l+1 + 1, . . . , m' r + 1 variables each with rrr! i polynomials of multidegree di for 
every 1 < i < r. We have rrr! x < Yljjti m 'j an d) for i 7^ 1, the condition rrii < Ylj^i m j 
implies that rrr! i < J2j=ii m 'jj therefore, renaming variables and polynomials, we are 
again under the assumptions of (III). 

5.2 A bound for the degree of the resultant 

This section is devoted to proving an upper bound for the degree D of the resultant of 
multi-linear polynomials appearing in our previous computations in terms of the total 
number n of strategies available to the players and the number 5 of totally mixed Nash 
equilibria of a generic game with the given structure. 

As before, we assume n = ri\ + ■ ■ ■ + n r with m, G N for every 1 < i < r. We 
consider the resultant TZ of a family of n + 1 multi-linear polynomials in r groups of 
ri\ + 1, . . . , n r + 1 variables each, consisting of a polynomial Fq of multi-degree do := 
(1, . . . , 1) and, for every 1 < i < r, a set of rij polynomials Fjp (1 < k < rii) of multi- 
degree di := (1, . . . , 1, 0, 1, . . . , 1), where the lies in the ith. coordinate. 
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Throughout this section, we will denote 



Si : 



= Bez 



(d , l;di,m; . ..; 



di,rii ~ 1 




i = l,...,r, 



the number of solutions to a generic polynomial system of the indicated multi-degrees. 
Recall that 5 = Bez nij ... )Tlr (di, rai; . . . ; d r ,n r ). 

(i) 

The resultant 7Z is a multihomogeneous polynomial in the coefficients of -Fo, Ft of 
degree 5 in the coefficients of Fq and Oj in the coefficients of for every 1 < i < r, 1 < 
k < rii (see, for instance, |25j). Therefore, the total degree of 1Z equals 



The following result provides an upper bound for 5i which will allow us to deduce an upper 
bound for D. 

Proposition 8 For every 1 < i < r , we have 5i < (n, + 1)5. 

Proof. Without loss of generality, we will prove the stated upper bound for i = 1. 

First, note that the formula in the right hand side of identity ([2]) for the compu- 
tation of the Bezout number of a multihomogeneous system of polynomial equations 
shows that this number is additive in each of the multi-degrees involved. Hence, as 5\ = 
Bez niv .. in ,.((fo, - l;d 2 ,n 2 ; . . . ;d r ,n r ) with d = d\ + ei, where e\ = (1,0,... ,0), 

we have 

5i = Bez ni) ... )nr (di,ni;d2,n 2 ; . . . ;d r ,n r ) + Bez nij ... jnr (ei, l;di,ni - l;d 2 ,n 2 ; . . . ;d r ,n r ) 
= 5 + Bez nir .. )flr (ei, l;di,ni - l;d 2 ,n 2 ; . . . ;d r ,n r ). 

Now, identity ([2]) implies that Bez„ 1) ... )n , r (ei, 1; d±, n\ — 1; d 2 , n 2 ; . . . ; d r , n r ) = where 
Zi = {(jn,---,jrn T ) / hi = l Jik iV(i,k) / (1,1) and #{jhk / jhk = = m V 1 < i < r}. 

In order to finish the proof, we will show that 4r~Zi < n i<5- Since 5 equals the cardinality 
of the set 3o introduced in we will compare the cardinalities of both sets 3i and 
3o- To this end, we define the following map from 3i to 3o : with a given n-tuple j := 
(1, ji2, ■ ■ ■ ,ji ni , ■ ■ ■ ,jrl, ■ ■ ■ ,jm T ) 6 Zl we associate the n-tuple j' G 5o which is obtained 
by exchanging the first coordinate of j (which equals 1) with the first one which is different 
from 1 and is located beyond the nith coordinate. Note that a necessary condition for two 
distinct n-tuples in 3i to lead to the same n-tuple in 3o by means of this assignment is 
that they coincide in all of their coordinates except for two of them located among the n\ 
coordinates n\ + 1, . . . , 2ni. Moreover, the vector consisting of these n\ coordinates must 
be of the form (1, . . . , l,jhk, • • • ) for both of them (possibly with no 1 at the beginning) 
and so, they can only differ in the length of the string of l's in this vector, which ranges 
between and n\ — 1. We conclude that each element of 3o is the image of at most rii 
elements of(l) 3i- It follows that #3i < ni#3o- 

Therefore, we have that 5± < 5 + m<5 = (ni + 1)5. □ 




(16) 



Ki<r 



26 



Using the previous result along with identity (fl~6j) for the degree D of the resultant, 
we conclude: 

Corollary 9 With the previous assumptions and notations, 



Remark 10 The bound stated in Corollary shows that all the algorithms presented in 
this paper are polynomial in the number of strategies n±, . . . ,n r of the r players, and the 
generic number 5 of totally mixed Nash equilibria of a game with the considered structure. 
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